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ABSTRACT 

The unusual mushroom-shaped HI cloud, GW 123.4-1.5, is hundreds of par- 
sees in size but does not show any correlations to HI shells or chimney structures. 
To investigate the origin and velocity structure of GW 123.4-1.5, we perform 
three-dimensional hydrodynamical simulations of the collision of a high-velocity 
cloud with the Galactic disk. We also perform a parameter study of the density, 
radius, and incident angle of the impact cloud. The numerical experiments in- 
dicate that we reproduce the mushroom-shaped structure which resembles GW 
123.4-1.5 in shape, size, position-velocity across the cap of the mushroom, and 
the density ratio between the mushroom and surrounding gas. GW 123.4-1.5 is 
expected to be formed by the almost head-on collision of a HVC with velocity 
~ 100 km s _1 and mass ~ 10 5 M about 5 x 10 7 yr ago. A mushroom-shaped 
structure like GW 123.4-1.5 must be infrequent on the Galactic plane, because 
the head-on collision which explains the mushroom structure seems rare for ob- 
served HVCs. HVC-disk collision explains not only the origin of the mushroom- 
shaped structure but also the formation of a variety of structures like shells, 
loops, and vertical structures in our Galaxy. 
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Introduction 



GW 123.4-1.5 revealed by the Canadian Galactic Plane Survey ( ITaylor et a/.ll2003l ) is 
a mushroom-shaped HI stru cture composed of a stem and a cap and is dissimilar to any 
common shells or c himneys (English et aZ.I l2000l ) . Using the kinetic distance of the Galac- 
tic rotation model, English et al\ (120001 ) derived properties of the mushroom cloud such as 
projected length, mass, velocity, mean excess column- density, and kinetic energy. The mush- 
room cloud extends a few hundred parsecs (~ 350 pc) at the assumed distance of 3. 8 ±1.2 kpc 
and has an HI mass of ~ 10 5 M Q , a mean excess column-density of Ah — 6 x 10 20 cm -2 , and 
a kinetic energy of ~ 2 x 10 50 ergs. As a model for the origin of the mushroom-shaped c loud 
GW 123.4-1.5, the gas bu oyancy model fenglish et aZ.lboOol : lAvillez k Mac Low IboOlh and 
the cloud collision model ( Kudoh fc Basu 1120041. hereafter KB2004) were proposed and nu- 
merical simulations were executed, respectively English et all ( 120001 ) carried out preliminary 
two-dimensional hydrodynamical simulations and showed that the origin of the mushroom- 
shaped GW 123.4-1.5 is the rise of buoyant gas from the explosion of a supernova above the 
Galactic midplane. Using three-dimensional hyd rodynamical s i mulations of large-scal e mod - 
eling of the interstellar gas in the galactic disk (lAvillez Il2000l ). lAvillez fc Mac Low I (120011 ) 
showed that the mushroom-shaped cloud can result from the buoyant rise of bubbles out of 
pools of hot gas created from a supernova or supernovae. However there is some difficulty in 
interpreting the origin of GW 123.4-1.5 by the buoyancy of a hot gas created by a single or 
multiple supernova explosions, because there are no observations of soft X-ray emission and 
IRAS sources near the base of GW 123.4-1.5. On the othe r hand, the collision of a h igh- 
velocity cloud (HVC) with the Galaxy, which was studied by iTenorio-Tagle et al\ ( 19861) for 



the fi rst time, is not restricted to the regions of active star formation (jTenorio-Tagle et al 
19871 ). KB2004 performed two-dimensional hydrodynamical simulations of the collision of 



a cloud with the Galactic disk and proposed that GW 123.4-1.5 is created by the impact 
of an intermediate- velocity cloud (IVCjj with the Galactic disk. Also they showed that the 
velocity structure across the cap of the mushroom-shaped structure in their simulation is 
consistent with that of GW 123.4-1.5. In the buoyant models, the velocity characteristic 
of the cap is not obvious. Unlike the gas buoyancy model, cloud-disk collision can easily 
explain the formation of GW 123.4-1.5 without inference for X-ray emissions. 

In this study we extend the two-dimensional hydrodynamical simulation of KB2004 
to three dimensions, in order to investigate the origin, dynamics, and velocity structure 
of a mushroom-shaped structure. Through numerical simulations and a parameter study, 



1 Historically, the clouds with \vi sr \ > 70 km s 1 are classified into HVCs and those with < 70 km s 1 are 
called intermediate- velocity clouds (IVCs). 
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we can reproduce the size, shape, and position-velocity of the mushroom-shaped structure 
which resemble those of GW 123.4-1.5. Consequently, we can infer (1) how GW 123.4- 
1.5 was formed in the Galactic disk and (2) why a mushroom-shaped structure is so rare 
on the Galactic plane. In the next section we describe our models and numerical method. 
Simulation results are presented in §111, followed by a summary and discussion in §IV. 



2. Numerical Simulations 



We simply suppose an adiabatic gas as in KB2004. As a galactic disk and halo model, 
we consider a hydrostatic equilibrium of the interstellar isothermal gas determined by the 
gravity and temperature. In this study, we use the Cartesian coordinates (x,y,z). The 
galactic disk is located at the centered x — y plane of a simulation box, and the z-axis is 
perpendicular to the galactic disk. Near the galactic midplane, density distribution of the 
interstellar isothermal gas under hydrostatic equilibrium is well approximated as 



p(z) = p exp [-(z/Hq) 2 ], 
where p is the density at the galactic midplane. The scale length is 

H = ^2/70^0, 

where 7 is the specific heat ratio 5/3 and c s q is the sound speed at the midplane. 
observed vertical gravitational field is app roximated as q 7 = —az, where a ~ 2/3 x 10" 



for th e solar neighborhood of the Galaxy fjSpitzer Ill978l ; iBahcall Ill984l ; iTenorio-Tagle et al 



I! 



(2) 



The 

29 s -2 



19871 ). We assume an external gravitational field and temperature distribution of interstellar 
gas at the initial stage as 

g z = -aH g t&nh(z/Hg), (3) 



T{z) = T + 0.5(T\ - T 



tanh 



\z\ - Zt 



(4) 



where H g is 1.5ifo an d T and Ti represent the temperatures of the disk and the halo. We 
assume T = 10 4 K, T\ = 10 5 K, z t = 0.5H Q , and za = 0.1H Q as those in the previous two- 
dimensional simulations (KB2004). As an initial condition for the Galactic disk and halo, we 
assume a hydrostatic equilibrium in the gas, which is determined by the Galactic gravity of 
equation (3) and the temperature distribution of equation (4). The density distribution of the 
hydrostatic equilibrium is nearly equal to equation (JT]) for z < z t , and it tends to be p(z) oc 
exp(—z/ Hi) for z > H g where H 1 = (Ti/T )('ycj l0 /aH g ). In addition to this hydrostatic 



equilibrium gas disk, we assume a cloud with density p c and velocity v c 
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Using three-dimensional velocity v c q = |v c | and the incident angle 8i, they are assumed as 
follows: 



Pail- 0.5 



1 + tanh 



I'd 



shift,- x w 



c0 



0.5 



v 



cosOj x t> 



c-0 



where r c = [(x — x c ) 2 + (y 



[z - z, 



yc - 

0.5 

\2lV2 



1 + tanh 

= 0, 

1 + tanh 



(5) 

(6) 
(7) 



is the distance from the cloud center (x c , y c , z c 



and v c0 = 100 km s 1 is a velocity of the cloud. v xc and v zc are parallel and perpendicular 



to the galactic disk, respectively. We take x c 



Vc 



0, Z r 



-2.5H Q , r wc = 0AH or 0.6H Q , 



rd = 0.1H , and r v 



+ ltd for a colliding cloud. 



We perform three-dimensional numerical simulation s using an Euleri an hydrodynamics 
code based on the total variation diminishing scheme (IRyu et q/.lll993l ). While radiative 
cooling and magnetic fields are likely to be important to the dynamics of interstellar gas in 
the galactic disk, we ignore them due to numerical complexity and limited computational 
resources. Radiative cooling and magnetic field will be considered in an upcoming paper. 
Simulations are made with 256 x 256 x 512 grid zones in x, y, and ^-directions respectively, 
and the computational box size is taken 0.7 kpc x 0.7 kpc x 1.4 kpc. A free boundary 
condition is imposed on all outer boundaries of the computational box. Simulations start at 
t — and last up to t cn d = 6io- All physical variables are expressed in units of the following 
normalization: Hq = 140 pc, c s q = 10 km s _1 , and to = Hq/c s q ~ 1.4 x 10 7 yr. 

This paper presents a total of seven simulations, which differ in density (p c ), radii (r c ), 
and incident angles of impact clouds. Model parameters are summarized in Table 1. 



Results 



3.1. Head-on Collision 



Firstly, we present the case of the head-on (9i = 0) collision to compare with the previous 
two-dimensional simulation (KB2004). Figure 1 and 2 show the evolution of the density and 



2 In eq. (5) of KB2004, '1 + tanh fedSiltZz 
same mistake is in eq. (6) of KB2004. 



was accidentally mistyped as 'tanh ( rc ^ R '*l rroc V. The 
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pressure as well as the velocity field on the y = plane in model HH, in which the same 
parameters are taken as those of KB2004; n = lcm -3 , r wc = 56(H /140 pc) pc. The top-left 
panel depicts the initial conditions. The cloud collides with the galactic disk at t/to ~ 0.2, 
and penetrates through the dense midplane at t/to ~ 0.3. At t/to — 0.8 we can see that 
a bow shock is traveling toward the halo ((x,z) ~ (0, 3.7ifo)) and a shell-like structure 
(density enhancement around (x,z) ~ (0,2.75i?o)) is formed by the cloud penetration (the 
third snapshot). Although the shell consists of a high-density gas, the pressure (see Fig. 2) 
is not especially high compared with the post-bow shock gas. This elongated shell structure 
originates in the tenuous halo gas and the dense disk gas. The inside of the elongated shell 
is pervaded by a gas coming from the impact cloud (t/t ~ 0.8 — 1.0). In the pressure 
plot (Fig. 2), two reverse shocks appear at t/to — 0.8. One is seen at the back of the cloud 
((x, z) ~ (0, 0)) and the other is around (x, z) ~ (0, 2H ) at the same time. Third and fourth 
snapshots of Figure 2 indicate both reverse shocks bear the signature of the Mach disk, which 
is formed by a reflection of the shock front on the 2-axis. On the other hand, in the galactic 
disk, a cavity is left after the cloud penetration. Three snapshots (from the fourth to the 
sixth) of Figure 1 indicate that the cavity is gradually being refilled with the galactic disk 
gas. This seems to be driven by the radially inward pressure force (t/to — 1-6 — 2.0). The 
reverse shock which is formed at t/to ~ 0.8 and propagates downward (seen at t/to ~ 1-6) 
interacts with the gas converging to the center of the cavity. The reverse shock hits the high- 
density gas at the age of t/to — 2.0, just when the cavity has been refilled. By the interaction 
between the shock wave and the converging gas, the converging dense gas is shock-heated 
and begins to rise along the z-axis (t/to ^ 2.2: 7th and 8th snapshots). The gas at high 
latitude is broken into eddies or vortices and forms the cap of a mushroom-shaped structure 
(t/to ^ 2.4). The shock-heated rising gas contributes to filling the central part of the stem 
which is seen in the 7th and the last snapshots (t/to — 2.4 — 2.8). Similarly to KB2004, 
the gas broken into eddies at the high latitude (cap) and the dense gas raised by the cavity 
replenishment (stem) have finally formed a mushroom-shaped structure. 

Figure 3 and 4 depict the evolution of the density and pressure as well as velocity 
field in model HL. In this model, the cloud is less dense no = 0.1cm -3 but larger r wc = 
84(if /140 pc) pc than that of model HH. Although the initial cloud masses assumed in 
models HH and HL are quite similar, the initial kinetic energy = (1/2) J pv 2 dV of the 
HVC in model HL is 40 % less energetic than that in the previous model HH. Because 
the density and velocity distributions of the cloud are not assumed to be uniform but to 
have smooth distributions (eqs. (5)-(8)). Dynamical evolution, density, and temperature 
distributions of the mushroom in model HL differ considerably from those in model HH. After 
a cloud collides with the galactic disk, a shell-like structure is formed inside the galactic disk 
at t/to — 1-0 (second and third snapshots). Even if the cloud masses in both models are 
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similar, the shell-like structure in model HL is different from model HH; the lateral diameter 
of the shell ~ 2.0H (third and fourth snapshots of Fig. 3) is much larger than that of the 
previous model ~ 1.0H (third snapshot of Fig. 1); in contrast, the vertical size is half that of 
model HH. In the third snapshot, we can see a reverse shock at (x, z) ~ (0, 1.25-ffo) moving 
toward the galactic plane at t/t ~ 1.2. This shock wave propagates and reaches the midplane 
at t/t ~ 1.8 (fourth snapshot). At the same time, the cavity is refilled with the disk gas 
approaching from lateral directions pushed by the interstellar pressure (t/t ~ 2.0 — 2.2). 
The wave interacts with the refilled dense gas and forms a high-pressure gas, which is seen 
near z ~ 0.5Ho in the 5th snapshot of Figure 4. This promotes the cavity replenishment and 
lifts the galactic dense gas along with the vortices (t/t > 2.2). In model HH, the vortex 
structures have been developed before the cavity replenishment has been finished. However, 
in model HL, the vortices are developed in the process of the cavity replenishment at the 
galactic plane. The vortices, of which the central location is rising along the stem, play a 
major role in forming the cap of the mushroom. Finally, the rising gas and vortices from 
the disk make a mushroom-shaped structure (t/to ~ 4.8 — 5.6). The rising gas is decelerated 
by the gravity of the galactic disk and eventually falls toward the disk, as seen in the last 
snapshot (t/t > 5.2). 

Figure 5 compares the density, velocity vector field (top panels), and temperature (bot- 
tom panels) of the mushroom-shaped structure on the y = plane in models HH (left) and 
HL (right). We can see the density and temperature differences between models HH and 
HL. In the case of model HH, the density shows a hollow in the stem which is filled with a 
hot gas. On the other hand, the stem of model HL has a high-density and low-temperature 
without any indication of the hollow. The difference comes from the formation processes 
of the stem. The stem of model HH was first formed from the gas which was originally a 
high-density disk and elongated into the low-density halo as the cloud collided with the disk 
and moved upward. Later on, a high-density and low-temperature gas blew up and occupied 
the inner hollow, after the disk was replenished and collided with a reverse shock. On the 
other hand, the stem of model HL is not formed until the disk gas blows up, because the 
disk gas is not sufficiently elongated upward when the cloud collides with the disk. Then, 
the stem of model HH has a hot hollow structure during its evolution, but the stem of model 
HL is always occupied by the dense and cold gas. This difference can also explain the reason 
why the mushroom structure of model HL appears later than that of model HH. 

KB2004 showed that the cloud density should be of the same order as the disk density 
to attain a mushroom-shape. Because the case of IVC (v zc = 50 km s" 1 ) showed a better 
agreement mainly with the height of the mushroom-shaped structure, they proposed that 
GW 123.4-1.5 was created by the impact of a high-density (no = 1 cm~ 3 ) IVC with the 
galactic gas disk. Our model HL shows that the mushroom-shaped structure can also be 
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made by the collision of a low-density (no = 0.1 cm -3 ) HVC (v zc = 100 km s _1 ) with the 
galactic disk, if the size of the cloud is sufficiently large. Whether for HVC or for IVC, the 
cloud-disk collision finally forms a mushroom-shaped structure. Therefore we may conclude 
that cloud collision with the Galactic disk is a promising model for the mushroom-shaped 
GW 123.4-1.5. 



3.2. Oblique Collision 

In this section, we explore oblique collisions with different incident angles (see Table 
1). As models of oblique collision, we assume the same cloud in model HL and change 
only the incident angle of the cloud. Figure 6 and 7 show the column-density and initial 
velocity (white arrows) of a cloud on the x — z plane and the y — z plane, respectively. We 
integrate the density along the y-axis (Fig. 6) and x-axis (Fig. 7) to get a column-density at 
t/t = 2.4 (middle panels) and t/t = 4.8 (bottom panels) for four models with different 
incident angles (#, is an angle between the vertical axis of the galactic disk and the initial 
velocity vector of a cloud), 5° (model 05), 10° (model O10), 20° (model O20), and 30° 
(model O30), respectively. The time evolution of model 03 (6*j = 3°) is shown in Figure 8 by 
the snapshot of column-density distribution. Figure 6 shows that a variety of structures such 
as mushroom-shape (bottom panel of model 05), shell/loop-like (middle panels of models 
O5-O30), and vertical structures (bottom panels of models 05-020) are effectively formed 
by cloud-disk collisions. After the cloud-disk collision, a variety of structures are quickly 
created by a subsequent process of cavity replenishment, gas rise, and eventual fall back 
due to the gravity of the galactic disk. The dynamical evolutions of oblique collisions are 
roughly similar except for the final structures. The resultant structures formed by cloud-disk 
collision are greatly affected by the incident angle. When we observe the cloud-disk collision 
on the y — z plane (Fig. 7), the projected structure (column-density integrating along the x- 
axis) evolutions in each model look very similar besides resultant structures (bottom panels). 
First, a shell and cavity are created, as seen at £/£ = 2.4 (middle panels), then change to 
vertical structures which penetrate the galactic disk vertically. Like the resultant structures 
observed on the x — z plane (bottom panels of Fig. 6), the vertical structures projected at 
x = +oo are mainly influenced by the incident angle at the x — z plane. From this parameter 
test, we can comprehend that cloud-disk collision of the large incident angle (model O10, 
O20, and O30) can hardly create a mushroom structure. Therefore if we can observe a 
mushroom structure along the x-direction, the incident angle at the x — z plane must be 
small 6i < 5°. 

Obliquely colliding clouds have initial velocity v xc and v zc given by equations (6) and 
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(8). The initial trajectory of the colliding clouds is represented by white arrows in the 
middle (t/t = 2.4) and bottom panels (t/t = 4.8) of Figure 6. At first glance, the resultant 
structures in the middle panels (especially for 0, > 10°) seem to be formed by the collision 
of a cloud moving from right to left, contrarily to the real situation. As a consequence of 
the cloud collision, a structure expanding in z < region is formed as well as the upwardly 
expanding main component. This countersplash, which is named for the structure seen in 
z < 0, seems to be a wake of the colliding cloud. If we misidentify this splash as a wake, we 
may consider the cloud moved from right to left. 

In a preliminary parameter test (we do not report the simulations in this paper), we 
considered several clouds which have a variety of velocity (~ 30 — 200 km s _1 ) and mass 
(~ 10 4 — 10 6 M Q ). When the kinetic energy of the cloud is larger than > 10 53 ergs, the 
structures after the cloud-disk collision are > 700 pc vertical structures. We did not explore 
these structures, because their size is much larger than the HI mushroom structure. On 
the other hand if the kinetic energy of the cloud is smaller than < 5 x 10 51 ergs {i.e., ~ 
30 km s -1 ), the colliding cloud does not penetrate the disk such as described in the previous 
two-dimensional simulation (KB2004) or does only make shell/loop-like structures near the 
galactic plane. In this study, our concern is focused on the formation of the mushroom- 
shaped cloud. We chose the incident kinetic energy of the cloud as ~ 10 52 ergs. Figure 6 
clearly shows that the cloud-disk collision can hardly make a mushroom-shaped structure, 
when the incident angle of a colliding cloud is larger than about 5°. Among our simulation 
results, the model 03 (6*; = 3°; Fig. 8) can describe GW 123.4-1.5 most closely. Then, the 
mushroom-shaped structure GW123.4-1.5 is likely produced by the approximately head-on 
collision of the cloud with a velocity of ~ 100 km s _1 , a mass of ~ 10 5 M Q , and thus a 
kinetic energy of ~ 10 52 ergs. 



3.3. A Model for GW 123.4-1.5 

Figure 8 shows the evolution of column-density on the x — z plane (integrating along 
the ?/-axis) for model 03. We can find the mushroom-shaped structures between t = 
56(t /14 Myr) Myr (at t/t = 4.0) and t = 67(t /U Myr) Myr (at t/t = 4.8). The 
height of the mushroom at the age of t = 67(to/14 Myr) Myr is about 350(i7o/140 pc) pc, 
which is the same size as that of GW 123.4-1.5. The lifetime of the mushroom structure is 
~ 10 7 yr and the width ratio of the cap to the stem is about 3:1. The mass of the cap is 
about twice that of the stem, and the total mass is estimated to be ~ 10 5 M Q . To obtain the 
mass of the mushroom-shaped cloud, we integrate the low-temperature gas with T < O.8X1. 
The mean column-density contrast between the mushroom and ambient medium is about 
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2:1. Although the mass ratio is a little smaller than the observed one 4:1 ([English et al 



20001 ) . most physical parameters of the mushroom-shaped structure in model 03 coincide 



with those of GW 123.4-1.5. 

In order to explain the velocity structure of the mushroom-shaped structure from the 
two-dimensional simulation, KB2004 assumed that the mushroom-shaped cloud is tilted to 
the line-of-sight and calculated the line-of-sight velocity from their two dimensional sim- 
ulation. However, without any presumption, we can reproduce the position-velocity map 
of the mushroom-shaped cloud from the three-dimensional simulation. Figure 9 shows the 
column-density and the position-velocity maps, which are equivalent to the observational 
longitude- velocity (/ — v) maps, for model 03. In order to distinguish the mushroom gas 
from the surroundings, we integrate only the low-temperature gas (T < O.8T1) along the 
y-axis (Fig. 9a) and the x-axis (Fig. 9b), respectively. Figure 9a corresponds to the case 
observed from y = —00 and Figure 9b from x = +00. The position-velocity plot of Figure 
9a indicates the symmetry with respect to v y = axis for both the bottom (z = 2AHq) 
and the upper (z = 3.2Hq) regions of the cap. On the other hand, the position-velocity 
plot of Figure 9b shows that the upper part of the cap has no significant velocity gradient, 
but the bottom has two distinct lobes blue-shifted with respect to the center of the cap 
by about 5 km s _1 . This velocity grad ient is very similar to the observed value of GW 



123.4-1.5 (Fig. lc of English et oZJl2000h . Observed (GW 123.4-1.5) and simulated (model 
03) mushroom-shaped clouds are compared in Table 2. Using model 03 which fits the 
observations best, we can infer that GW 123.4-1.5 is created by a HVC with a velocity of 
~ 100 km s™ 1 coming toward us with vi sr ~ 5 km s _1 and a mass of ~ 10 5 M collided with 
the Galactic disk about 5 x 10 7 yr ago. 



4. Summary and Discussion 

We perform three-dimensional hydrodynamical simulations for the impact of a HVC 
with the Galactic disk, in order to explore the formation of the mushroom-shaped HI cloud 
GW 123.4-1.5. The main results can be summarized as follows: 

• A mushroom-shaped structure can be formed not only by a high-density IVC collision 
with the Galactic disk (KB2004) but also by a low-density HVC collision. 

• The resultant structures formed by cloud-disk collision are greatly affected by the 
density and incident angle of a colliding cloud. 

• GW 123.4-1.5 is expected to be formed by the almost head-on collision of a HVC with 
velocity ~ 100 km s _1 and mass ~ 10 5 M & about 5 x 10 7 yr ago. 
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Recently, lAsgekar et all (120051 ) have tried to find a mushroom-shaped HI structure 
like GW 123.4-1.5 in the Canadian Galactic Plane Survey. However, they did not find 
other mushroom shape structures. Our simulations could explain why mushroom-shaped 
structures are so infrequent in our Galaxy. A nearly head-on collision is needed to form a 
mushroom-shaped structure. In a rotating disk, the head-on collision of a cloud with the 
Galactic plane would be rare in a general way. 

HVC-disk collision can not only explain the origin of a mushroom-shaped structure 
but also the formation of a variety of structures like shells, loops and vertical structures 
in our Galaxy. For example, in model O30 (see the middle panel of Fig. 6), we can find 
an interesting structure, a loop which resembles that found by NANT EN in the central 
molecular zone within ~ 1 kpc from the Galactic center (IFukui et a/.ll2006l ). The two loops of 
CO emission obtained with NANTEN have strong velocity gradients in a longitude- velocity 
di agram: ~ 80 km s" 1 per 250 pc and ~ 60 km s _1 per 150 pc in each loop (see Fig. 2 



of IFukui et al\\200w . Although the projected velocity distribution of our simulation (~ 
20 km s" 1 ) is somewhat smaller than that of observation (~ 60 — 80 km s _1 ), the size of the 
loop, ~ 300(ifo/140 pc) pc, is si milar to that o f the molecular loop. In order to explain the 
formation of the molecular loop, Fukui et al. J2OO6}) pr oposed a magnetic flotation model 



Parker 



(119661 ) instability). From our preliminary 



(due to magnetic buoyancy caused by the 
result of model O30, we also propose that HVC-disk collision is another possible model of 
formation of the molecular loop structure at the Galactic center. This issue will be studied 
in detail in an upcoming paper. 

The present study as well as the previous study (KB2004) are restricted t o hydrodynam 



ical s i mulations. When the magnetic fie ld parallel to the Galactic disk exists (jSantillan et al. 



20001 ; iFranco et a/.ll2002l ; iNozawa 1120051 ) , we also expect the magnetic buoyant effect (Parker 



instability) . Both the impact of the cloud and the Parker instability (jSantillan et q/.lll999 



20041 . 120071 ) may create complicate structures in the Galactic disk such as worms, loops, 
shells, and chimneys. Moreover the interaction between the infalling HVCs and the magne- 
tized disk must be an elementary p rocess in the interstell ar dynamics. Using high-resolution 
two-dimensional MHD simulation, ISantillan et all ( 120071 ) pointed out a possibility that the 
rain of compact HVCs onto the disk can maintain transonic turbulent motion in the warm 
phase of HI. In this sense, we believe that high-resolution three-dimensional MHD simula- 
tion will broaden our understanding of the structures formed by HVC collisions and bring 
us more interesting phenomena to be important in the interstellar dynamics. In the near 
future, we will perform three-dimensional magnetohydrodynamical simulations. 
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Table 1. Model Parameters for Simulations 



Model 


n cQ ( cm 3 ) a 




^(deg) b 


HH 


1.0 


0AH 


0.0 


HL 


0.1 


0.6H 


0.0 


03 


0.1 


0.6H 


3.0 


05 


0.1 


0.6H 


5.0 


010 


0.1 


0.6H 


10.0 


020 


0.1 


0.6H 


20.0 


030 


0.1 


0.6H 


30.0 



a n c0 is the number density at the center 
of the HVC. 

h 9i is an angle between the vertical axis 
of the galactic disk and the initial vector 
of a cloud. 
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Table 2. Comparison of observed (GW 123.4-1.5) and simulated (model 03) 

mushroom-shaped cloud 



GW 123.4-1.5 £ 



model 03 



height ( pc) 
total mass (10 5 M Q ) 
mass ratio (cap:stem) 
Av c ( km s _1 ) 

mean column-density contrast (mushroom: ambient medium) 



350 b 
1.55 
3.4:1.0 

5 

4:1 



350 • (#o/140 pc) 
0.76 
1.7:1.0 

5 

2:1 



'Observational data from English et all (120001 ) . 



b A distance of 3.8 kpc was assumed in English et al\ ( 120001 ). 
c Velocity difference between the lobes of the cap and the central cap. 



Fig. 1. — Time sequence of density distributions (in logarithmic scale) and velocity field at 
the y = plane at t/t = 0,0.2,0.8, 1.2 (top:left to right), 1.6,2.0,2.4, and 2.8 (bottomdeft 
to right) for model HH. The abscissa represents the x-axis and the ordinate represents the 
z-axis. The longest arrow correspond to about 100 km s _1 . 




Fig. 2. — Time sequence of pressure at the y = plane at t/t 
right), 1.6,2.0,2.4, and 2.8 (bottom:left to right) for model HH, 



= 0,0.2,0.8,1.2 (topdeft to 




Fig. 3. — Same as Fig. 1, but for model HL. Time sequence is t/to = 0,0.4, 1.2, 1.8 (top: left 
to right), 2.2,3.2,4.8, and 5.6 (bottomdeft to right). 




Fig. 4.- 



Same as Fig. 2, but for model HL. 




Fig. 5. — Density, velocity (top panels) and temperature distributions (bottom panels) on 
the x — z plane, (a) is for model HH at t/to = 2.4 and (b) is for model HL at t/to = 4.8. 
Temperature (in logarithmic scale) is normalized in 10 4 K. 
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Fig. 6. — Column-density on the x — z plane (integrating along the y-axis) at three epochs 
(t/to = (top panels), 2.4 (middle panels), and 4.8 (bottom panels), respectively), for four 
models 05 (0< = 5°), O10 (10°), O20 (20°), and O30 (30°). At t/t = (top panels), initial 
velocity fields on the y = plane are shown. White arrows in the middle and bottom panels 
represent the initial trajectory of the cloud. The image scales are logarithmic for o = J pdl. 
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Fig. 7. — Same as Fig. 6, but on the y — z plane (integrating along the x-axis). 
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Fig. 8. — Time evolution of column- density of cool gas on the x — z plane. The snapshots are 
made at t/t = 0,0.8,1.6,2.4 (top:left to right), 3.2,4.0,4.8, and 5.6 (bottomdeft to right) 
for model 03. 
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Fig. 9. — Column-density (left panels) and position- velocity (right panels) at t/t = 4.8 in 
model 03. Position-velocity is measured for the lower (z = 2AH ) and upper (z = 3.2H ) 
regions of the caps on (a) the x — z plane (integrating along the y-ecxis) and on (b) the y — z 
plane (integrating along the x-axis). 



